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Abstract 

A one-dimensional model of a gel-fueled rocket com- 
bustion chamber has been developed. This model includes 
the processes of liquid hydrocarbon burnout, secondary 
atomization, aluminum ignition, and aluminum combus- 
tion. Also included is a model of radiative heat transfer 
from the solid combustion products to the chamber walls 
Calculations indicate that only modest secondary atomiza- 
tion is required to significantly reduce propellant burnout 
distances, aluminum oxide residual size and radiation heat 
losses. Radiation losses equal to approximately 2-13 % of 
the energy released during combustion were estimated. A 
two-dimensional, two-phase nozzle code was employed to 
estimate radiation and nozzle two-phase flow effects on 
overall engine performance. Radiation losses yielded a 1 
% decrease in engine Results also indicate that sec- 
ondary atomization may have less effect on two-phase 
losses than it does on propellant burnout distance and no 
effect if oxide particle coagulation and shear induced drop- 
let breakup govern oxide particle size. Engine 1^ was 
found to decrease from 337.4 to 293.7 seconds as gel 
aluminum mass loading was varied from 0-70 wt%. En- 
gine Igp efficiencies, accounting for radiation and two- 
phase flow effects, on the order of 0.946 were calculated 
for a 60 wt% gel, assuming a fragmentation ratio of 5. 

Nomenclature 

a- absorption coefficient 

A- area 

c p - specific heat 

C D - (frag coefficient 

d- diameter 

F n (T)- fraction of total radiation emitted by a blackbody 

at temperature T in the wavelength interval n 
h- average convective heat transfer coefficient 
i- specific enthalpy 

1- radiant intensity 

I 0 ,l|- modified Bessel functions 

j- droplet size class 

k- thermal conductivity 

L- number of wavelength intervals 

m- mass 

rh- droplet mass flux or mass flux through the com- 
bustion chamber 

M- number of droplet size classes 

n- wavelength interval 
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Nu- Nusselt number 

N/t- number of droplets passing through combustor 

control volume per unit time 
r s - droplet/particie surface radius 

r- radial coordinate/position in terms of optical 

thickness 

R- chamber radius 

R- chamber radius in terms of optical thickness 

q r - radiation heat flux from condensed combustion 
products to wall 
t- time 

T- temperature 

Up- particle velocity 

x- axial location in the combustor model 

Greek 

a t ,a 2 - control switches 

P* fragmentation ratio (number of secondary droplets 
produced per initial droplet) 

<j>- angular coordinate 

r|- fraction of aluminum droplet mass that forms as 
A1 2 0 3 residual on droplet surface 
X- wavelength 

p- density 

o- Stefan-Boltzmann constant 

a s - scattering coefficient 

0- angular coordinate 

Q 0 - scattering albedo 

Subscripts 

Al- aluminum 

b- boiling, propellant burnout or blackbody 
f- aluminum oxide fume particles or aluminum 

flame 

fg- liquid-gas 

g- gas 

LH- liquid hydrocarbon 

o- initial 

Ox- Aluminum oxide residual 

p- particle/droplet 

w- wall 

Introduction 

Gel propellants, in which fine particles are suspended 
in a gelled combustible liquid, offer potential performance 
and/or safety advantages over conventional liquid and 



solid propellants in rocket applications Theoretical per- 
formance evaluations show that gel propellants may pro- 
vide increases in specific impulse <l sp ) and/or propellant 
density over conventional liquid propellants 1-7 A number 
of potential Uqutd/solid constituent combinations were 
theoretically evaluated over a range of mixture ratios and 
solid constituent loadings for use in a LEO-GEO mission 
with fixed propellant tankage volume and vehicle dry -mass 
with no limit on initial gross mass 7 Based on propellant 
toxicity and vehicle payload considerations, aluminum was 
predicted to be the best candidate solid constituent. 

More detailed analyses of earth-to-LEO, LEO-to-GEO. 
LEO-to-Moon. LEO-to-Mars. and several planetary' fly-by 
missions, incorporating appropriate total vehicle 
mass/volume limitations and propellant density effects on 
vehicle dry -mass, also were performed to compare various 
aluminum gels with their liquid bi-propellant counter- 
parts. 3-6 These studies indicate that gel propellants are 
likely to provide the greatest benefit on earth-to-LEO. and 
high-energy planetary missions with NTO/MMH/A1 gels 
also providing significant performance gains over 
NTO/MMH in LEO-to-GEO and LEO-to-Moon missions. 
For example, replacing the space shuttle solid boosters 
with an RJM/O 2 /AI or NTO/MMH/ A1 propellant combi- 
nation. while maintaining the same booster dimensions, 
could theoretically increase maximum payload by 14-35% 
over the current shuttle payload. 6 Similarly, the addition 
of aluminum to nitrogen tetraoxide/ monomethyl hydra- 
zine (NTO/MMH) may permit planetary missions that 
would be impossible with neat NTO/MMH. 5 

In the area of propellant rheology, recent gains have 
been made in the development of stable gels that are easily 
sprayed. 8 ' 10 but more work in this area is still required. In 
related work, testing has demonstrated that expulsion 
pumping is feasible in a small engine system and that gel- 
fueled engines can be throttled over a range of operating 
conditions, shut-down, and restarted as desired. 11 In ad- 
dition. experimental results indicate that using aluminum 
can reduce combustion instabilities associated with liquid 
hydrocarbon propellants. 2 

Although the above investigations are promising, it 
should be noted that the theoretical performance studies 1-7 
do not incorporate 1^ losses associated with aluminum 
combustion such as increased propellant combustion times, 
radiation heat transfer from condensed combustion prod- 
ucts to the chamber walls, and nozzle two-phase flow 
losses due to these same condensed products. Because I S p 
efficiency losses of only 1.5^% may be sufficient to elimi- 
nate the payload benefits of using aluminum gels, minimi- 
zation of aluminum combustion related losses is desirable. 
In addition, these performance losses must be determined 
before the performance of gel propellants can be accurately 
evaluated. 

Since previous work 1213 has shown that the individual 
aluminum particles in a gel droplet can form an agglomer- 
ate that bums as a single aluminum droplet it is apparent 
that small gel droplets are beneficial in reducing perform- 
ance losses. Fine atomization of gel propellants, however, 


can be difficult to achieve due to their v iscous nature For- 
tunately. secondary atomization. 14-19 in which a droplet 
spontaneously shatters into a number of smaller droplets 
due to internal vaporization of the liquid earner, may pro- 
duce the desired small droplets 

Although secondary' atomization has been proposed as 
a means of reducing gel-related performance losses, previ- 
ous research has focused on the secondary atomization 
mechanism itself rather than its effects on engine perform- 
ance. The focus of our present inv estigation, therefore, is 
to theoretically examine the effects of secondary atomiza- 
tion on engine performance losses and. incorporating these 
performance loss estimates, evaluate the performance po- 
tential of gel-fueled engines. To accomplish these objec- 
tives, a one-dimensional combustor model has been de- 
rived to evaluate secondary atomization effects on propel- 
lant burnout distances and radiation heat losses in a rocket 
combustion chamber. This combustor code is used in 
conjunction with the Solid Propellant Rocket Motor Per- 
formance Prediction Computer Program (SPP), 20 a two- 
dimensional. two-phase rocket code, to provide a prelimi- 
nary evaluation of nozzle two-phase flow losses and 
overall engine performance. 

Combustion Chamber Model 
One-Dimensional Model Description 

In brief, a radially uniform spray, consisting of four 
droplet size classes, enters the combustion chamber and 
burns in a process incorporating liquid carrier burnout, 
droplet secondary atomization, aluminum agglomerate 
heat-up and combustion, two-phase particle flow, and ra- 
diation heat transfer from solid combustion products to the 
chamber walls. The propellants used in the combustor 
model are a JP-10/A1 gel and a preheated gaseous 0 2 oxi- 
dizer. JP-10, a pure hydrocarbon (C 10 H 16 ), was chosen as 
the gel hydrocarbon component instead of RP- 1 to avoid 
the complexity of modeling multi -component droplet com- 
bustion. The combustor flow is modeled using a single 
product-phase containing both gases and small A1 2 0 3 
fume particles, and three additional flow phases for each 
droplet size class: a liquid hydrocarbon phase, an alumi- 
num phase, and a phase containing large A1 2 0 3 residuals 
The A1 2 0 3 must be separated into two flow phases because 
of the two oxidation mechanisms inherent in aluminum 
combustion. 21 ’ 22 In the first oxidation mechanism, large 
A1 2 0 3 residuals are formed through droplet surface con- 
densation/oxidation, and in the second, very small A1 2 0 3 
fume particles are produced through vapor-phase oxida- 
tion. Since the A1 2 0 3 residuals are much larger than the 
fume particles and are attached to the surface of the alu- 
minum droplets it is necessary to model the A1 2 0 3 residu- 
als separately from the fume particles. Including the fume 
particles in the gas-phase flow requires the assumptions of 
no temperature or velocity slip between the fume particles 
and the gases, which are reasonably valid, and greatly 
simplifies the problem solution. 

As an A1 droplet bums, it continually produces A1 2 0 3 
fume particles at the oxide boiling point which subse- 
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qucntly equilibrate with the gas-phase temperature and 
velocity. This process results in a range of fume panicle 
temperatures and velocities as new panicles are created 
and older particles continue to equilibrate with the gases. 
If the temperature and velocity slip were not neglected, 
many additional flow phases would be required to accu- 
rately model this range of fume particle temperatures and 
velocities. 

Mass Conserv ation: Using the definition of velocity. 
u=dx/dt. to relate the differential variables, dx and dt. the 
following equation can be wntten for system mass conser- 
vation for M droplet size classes: 
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dm dXH dm dA1 | dm d0x 'j 
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where it is assumed that no mass is added to or removed 
from the chamber except at the injector face and the 
chamber exit, and that steady -state operation prevails. 

The three time derivatives on the right hand side of 
Eq. L are found from hy drocarbon droplet gasification and 
aluminum droplet combustion models which are discussed 
later in this paper. The value of N/tj for each droplet size 
class can be determined from the total initial gel mass flux 
and a normalized droplet size distribution. 

Energy Conservation: The steady-state energy bal- 
ance for the system, including radiation heat losses, can be 
expressed as 
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Since gas flow optical properties are dominated by the 
small A1 2 0 3 fume particles . 23 the radiation term is inde- 
pendent of the individual droplet size classes and is not 
included in the summation term. Each of the terms on the 
right hand side of Eq. 2 can be represented as follows by 
expanding the derivatives and assuming that the only 
changes in specific enthalpy are due to convective heat 
transfer with the gaseous product-phase: 
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The coefficients and a^, equal to 0 or 1, are used to 
provide the correct terms, depending on whether the liquid 
hydrocarbon, aluminum, or A1 2 0 3 mass flux is being ex- 
amined. In the hydrocarbon vaporization model, the bulk 
droplet temperature is assumed to remain at the initial 
droplet temperature while only a thin surface layer of hy- 


drocarbon is heated to the hydrocarbon boding tempera- 
ture before being vapon/ed Since the liquid hydrocarbon 
enthalpy does not vary with chamber axial location. 0^=0 
and a 2 =l. In this case, the enthalpy, i. equals the initial 
liquid hydrocarbon specific enthalpy, and dn^/dt is the hy- 
drocarbon mass vaporization rate from a single droplet. 

Following hydrocarbon burnout, an agglomerate of 
aluminum particles remains 12 ,3 The aluminum agglom- 
erate temperature rises from the hy drocarbon boiling tem- 
perature. through the aluminum melting point, to the 
aluminum ignition temperature as heat is transferred from 
the gas flow r to the agglomerate, at which point aluminum 
combustion commences. In the combustor model, this 
agglomerate ignition/heat-up process is approximated as 
convective heat transfer from the gas flow to a spherical, 
uniform-temperature agglomerate. Heat transfer to the 
agglomerate causes the agglomerate temperature to rise 
from the hydrocarbon boiling temperature to the aluminum 
melting temperature. The agglomerate temperature is then 
held constant until sufficient energy' to melt the entire ag- 
glomerate has been transferred from the gas-flow, at which 
point it is assumed that aluminum combustion begins. 
Heat-up of the molten aluminum droplet to the aluminum 
boiling temperature is accounted for in the aluminum 
combustion model, as discussed below. Since it is as- 
sumed that no aluminum vaporization occurs during ag- 
glomerate heat-up and melting, and because the specific 
enthalpy does change, a^l and 0^=0 during these stages 
of the aluminum ignition/ combustion process. Once the 
agglomerate melts, however, the aluminum specific 
enthalpy remains constant for the same reasons as dis- 
cussed for the liquid hydrocarbon, and combustion begins 
to occur. Chiring this portion of the ignition/combustion 
process, otj=0, 0 ^= 1 , and the specific enthalpy, i, equals 
that of liquid aluminum at its melting temperature 
Since the A1 2 0 3 residual forms on the aluminum drop- 
let surface, the A1 2 0 3 residual is maintained at the alumi- 
num boiling temperature, and therefore maintains a con- 
stant specific enthalpy, as long as any aluminum remains 
in the droplet. For the time period prior to aluminum 
burnout, 04 = 0 , 04 = 1 , and the specific enthalpy is that of 
A1 2 0 3 at the aluminum boiling temperature. Once alumi- 
num combustion and A1 2 0 3 formation ceases. a { =i. a 2 = 0 , 
and the particle temperature is allowed to equilibrate with 
the gas temperature. The total system energy balance is 
represented by the substitution of Eq. 3 with appropriate 
values of c^, a^, and i into Eq. 2 

Radiation Heat Transfer: Radiation from the solid com- 
bustion products to the chamber wails is a participating 
media phenomenon, requiring a solution of the radiative 
transfer equation . 24 Expressed in cylindrical coordinates 
this equation is 
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where I=I(r,B,4>). It should be noted that the radial posi- 
tion. r. is based on optical thickness and is not a physical 

position (i.e., r = J(a + a s )dr) The scattering albedo, 
is defined as 


where a and a s are determined from Mie theory. In the 
Mie scattering calculations the fume particles are assumed 
to be uniform in size ( 1 nm), and the refractive indices are 
those of A1 2 0 3 smoke. 23 

Since the flow is likely to be optically thick, 23 25 a cy- 
lindrical P { diffusion approximation is used to simplify the 
radiative transfer equation. By assuming that I b f is inde- 
pendent of radial location, the following equation for ra- 
diative heat flux is obtained: 25 



(6) 


The variable ^ is based on the flow optical properties 
and is defined as 
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If ^(k) is assumed to be constant over a wavelength in- 
terval. Eq. 6 can be expressed as the following summation: 
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Momentum Conservation: The gas-phase momentum 
equation in the combustion chamber is trivial assuming a 
negligible chamber pressure gradient and no body forces. 
Neglecting virtual mass and Bassett forces, the particle 
momentum equation for a given droplet size class can be 
expressed as 


du p _2 ~ u Jj u , ~ u ,. 

d * 4 PpdpUp 


Gel Combustion and Secondary Atomization: The 
physics of gel droplet combustion and probable secondary 
atomization processes are described in other work. 12-19 In 
the present model, rigid-shell induced secondary' atomiza- 
tion 17-19 is assumed. In the combustor code, secondary' 
atomization is presumed to occur when the droplet diame- 
ter reaches the predicted ngid-shell diameter. Although 
secondary atomization actually occurs some time after 
rigid shell formation, the time interval is currently un- 
known and is therefore neglected in the combustor model 
Particle size distribution after secondary atomization is 
also presently unknown and is treated as a system variable. 
In the combustor model, a droplet undergoing secondary 
atomization is assumed to shatter into a specified number 
of equal-size secondary droplets Defining the fragmenta- 
tion ratio. (3, as the number of secondary droplets produced 
per initial droplet, a new value of N/ij can be expressed as 



(10) 


Equation 10 is then used to determine a new droplet di- 
ameter, dj, given the gel mass flax of the size class 
Hydrocarbon Vaporization. Based on the relative prox- 
imity of the gel droplets to each other, it is unlikely that 
the droplets are surrounded by individual flames There- 
fore, hydrocarbon combustion is represented using a 
spherically symmetric droplet evaporation model with 
convective effects incorporated through film theory 26 As 
mentioned previously, droplet heat up is approximated by 
assuming that heat transfer from the gas flow only affects a 
thin liquid layer at the droplet surface. Assuming quasi- 
steady, spherically symmetric droplet vaporization with 
unity Lewis number, uniform droplet temperature, and 
constant thermophysical properties, the hydrocarbon va- 
porization rate can be expressed as 


d m d.LH 
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Aluminum Combustion. As mentioned previously, 
aluminum combustion proceeds through two oxidation 
mechanisms: aluminum vapor oxidation and droplet sur- 
face condensation/oxidation. This aluminum combustion 
process is treated using a simple combustion model in 
which a specified fraction of the vaporizing aluminum is 
involved in the surface oxidation/condensation mecha- 
nism. It is assumed that all heat released by this surface 
mechanism goes into the droplet with none lost to the sur- 
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rounding gases Using the same assumptions as the hydro- 
carbon model and assuming that A1 2 0 3 does not interfere 
with aluminum vaporization, the droplet vaporization rate 
can be expressed as 


dt 


27ik ? NurJ 


c pj(Tf T'b.Ai) 


J T c p.AidT-ntfg.ox 


1 fg.Al + 


( 12 ) 


It should be noted that i fg is not a true enthalpy of 
vaporization since A1 2 0 3 does not exist in a vapor state. 
Therefore. i ft>Al 0j {S actuary the enthalpy change for the 

chemical reaction, 2Al (g) +3/20 2 ^Al 2 0 3(l y occurring at 
the droplet surface temperature. Droplet lifetime compari- 
sons with an empirical correlation for aluminum droplet 
combustion in solid rockets 20 demonstrated good agree- 
ment. Average thermophysical properties are estimated as 
recommended in previous droplet combustion work. 27 

Solution Method : Equations 1,2,9, 1 1, and 12 form the 
governing equation set and are numerically integrated in 
the axial direction using the IMSL Dverk integration rou- 
tine. 28 Gas-phase temperature, density, and composition 
are calculated using the ST AN JAN chemical equilibrium 
subroutine. 29 The temperature dependent gas-phase con- 
ductivity and viscosity are approximated as those of 0 2 to 
avoid the complexities of calculating gas mixture proper- 
ties. Details of these and other properties used in the com- 
bustor model can be found in reference 30. 

One-Dimensional Combustor Model Results 

The one-dimensional combustor code was exercised us- 
ing the conditions presented in Table !. These values were 
chosen to simulate an upper-stage booster and represent 
the maximum 1^ operating point for a 60 wt% aluminum 
gel. Because propellant mass flow rate in a rocket is gov- 
erned in part by nozzle geometry, the SPP nozzle code 20 
was used to determine the propellant mass flow rates. 

Since an appropriate spray distribution remains to be de- 
termined from gel atomization research, the arbitrary nor- 
malized droplet size distribution presented in Fig. 1 was 
used. 

Gas temperature and composition, assuming a frag- 


Table 1. 

Model Operating Conditions 

Chamber Diameter 

0.19 m 

Chamber Pressure 

38 atm 

Liquid Carrier 

JP-10 

Gel Aluminum Loading 

60% 

Gel Flow Rate 

10.78 kg/s 

Oxidizer Flow Rate 

11.75 kg/s 


O 

z 



DROPLET SIZE (yun) 


Figure 1. Normalized droplet size distribution used in the 
combustor code. 


mentation ratio, p, of 5 are presented in Fig. 2. The jag- 
ged shapes of both the temperature and composition pro- 
files arises from using only four droplet size classes instead 
of a continuous droplet size distribution. The slow rate of 
temperature increase in the region between 0 1 and 0 25 m 
is caused by the large enthalpy transfer from the gas flow 
to heat the aluminum after hydrocarbon burnout 

In Fig. 3, propellant burnout distance is plotted versus 
fragmentation ratio, p, to illustrate the potential benefits of 



AXIAL LOCATION (m) 

Figure 2. Gas temperature and major species mass fluxes 
versus axial location. Data are for a 60 wt% aluminum gel 
assuming a fragmentation ratio of 5. 
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Figure 3. Propellant burnout distance as a function of sec- 
ondary atomization fragmentation ratio for a 60 wt% 
aluminum gel. 


Figure 4. Final oxide residual diameter as a function of 
secondary atomization fragmentation ratio for a 60 wt% 
aluminum gel 


secondary atomization in reducing required engine resi- 
dence times. It is readily apparent that only slight secon- 
dary atomization is required to significantly reduce propel- 
lant burnout distance, and that higher secondary atomiza- 
tion intensities, represented by larger fragmentation ratios, 
have a lesser marginal effect, primarily due to the fact that 
droplet lifetime is inversely proportional to droplet surface 
area, which increases as (p) 2/3 . 

Final oxide residual diameter as a function of fragmen- 
tation ratio is shown in Fig. 4. Similar to the trend seen in 
propellant burnout, small secondary atomization intensi- 
ties significantly reduce final A1 2 0 3 residual diameter, 
with greater atomization intensities providing decreasing 
marginal reductions in residual diameter. Final residual 
particle diameter, however, is proportional to (l/p) 1/3 in- 
stead of (1/P) 2 3 as is droplet lifetime. 

Radiation losses from condensed combustion products 
to the chamber walls as functions of fragmentation ratio 
and aluminum mass loading are shown in Fig. 5. Cham- 
ber operating conditions are those that produce maximum 
1^ for each aluminum loading. Radiation losses were 
found to be strongly influenced by the fragmentation ratio, 
primarily due to changes in propellant burnout distance. 
Increasing the gel aluminum mass loading for a given 
fragmentation ratio was also found to increase radiation 
losses because of increased chamber temperature and 
slight increases in propellant burnout distance with in- 
creasing aluminum mass loading. Radiation losses were 
approximately equal to 2-13 % of the energy released 
during combustion, depending on the aluminum mass 
loading and fragmentation ratio. 

It should be cautioned that the above radiation heat 
transfer calculations are only approximate. Axial radia- 


tion heat transfer, which has been neglected, may alter the 
droplet combustion process and/or the actual radiation 
losses. In addition, research 29 has indicated that the opti- 
cal properties of condensed combustion products can vary 
considerably depending on propellant composition and 
oxidizer/fuel mixture ratio. 



Figure 5. Radiation heat loss as a percentage of combus- 
tion energy released as a function of fragmentation ratio 
and gel aluminum mass loading. 
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O/F MIXTURE RATIO 

Figure 6. Engine I sp as a function of oxide particle size 
model and O/F mixture ratio for a 60 wt% aluminum gel. 
Dashed line indicates inclusion of radiation losses. Model 
assumes an Extended Delta nozzle profile. 


SPP Modeling of Engine Performance 
Nozzle Two-Phase Flow 

Two-phase flow effects in the engine nozzle were esti- 
mated using the SPP rocket code. 20 Nozzle inlet condi- 
tions were determined using a chemical equilibrium mod- 
ule of the SPP code that provides a mechanism to account 
for the radiation losses predicted by the one-dimensional 
combustor code. 

Because particle interaction and mass transfer between 
the particles and gas are neglected in the SPP code, the 
A1 2 0 3 residual size distribution predicted by the combustor 
code may be incorrect for nozzle performance calculations. 
Since the A1 2 0 3 particles are molten throughout most of 
the nozzle, and because small particles accelerate more 
quickly than large particles, particle size may increase 
through coagulation. Similarly, additional A1 2 0 3 may be 
produced through the recombination of gas-phase radicals 
as the exhaust gases cool during expansion, resulting in 
the nucleation of additional particles and/or growth of 
previously formed particles. Particle size may also de- 
crease due to shear breakup of droplets, particularly in the 
throat region of the nozzle. 

Because of these uncertainties in A1 2 0 3 particle size, 
two methods of estimating particle size, which should 
bound the true particle size, are used in the evaluation of 
nozzle performance. In the first method, we use an A1 2 0 3 
particle size distribution that is determined by the one-di- 
mensional combustor code, making secondary atomization 
the primary mechanism governing particle size. In the 
second method, secondary atomization is assumed to have 
no effect on mean particle size; rather, coagulation, parti- 



O/F MIXTURE RATIO 

Figure 7. Engine 1^, as a function of O/F mixture ratio 
and gel aluminum mass loading. Dashed lines represent 
inclusion of radiation losses. 


cle surface growth, and shear induced droplet breakup are 
assumed to be the dominant mechanisms affecting particle 
size. Assuming that these coagulation and breakup 
mechanisms are comparable to those in solid nozzles, the 
following solid motor correlation of mass median A1 2 0 3 

particle diameter, as a function of nozzle throat di- 
ameter was employed; 20 

D* = 3.63D, 2932 . (13) 

This results in a particle mass median diameter of 5.6 pm. 
Engine Performance Results 

A comparison of two-phase flow effects on 1^ for the 
above methods of determining A1 2 0 3 particle size are pre- 
sented in Fig. 6 for a 60 wt% aluminum gel. As seen here, 
secondary atomization may reduce two-phase flow losses 
but not as significantly as propellant burnout distance (4% 
versus 62% decrease, cf. Fig. 3). To illustrate the separate 
contributions of radiation and two-phase flow losses, ra- 
diation losses were incorporated in the case represented by 
the dashed line. Given a fragmentation ratio of 5. radia- 
tion losses yield a decrease in L, of approximately 1 % 
compared to the 4 % resulting from two-phase losses. 

Using the solid motor correlation for particle size. 1^, 
was calculated fin a range of aluminum mass loadings and 
propellant mixture ratios and compared with 1^ calcula- 
tions for a JP-10/O 2 bi-propellant. Figure 7 shows that 1^ 
decreases with increasing aluminum loading and that the 
maximum 1^ mixture ratio becomes richer, as has been 
predicted by other studies. 2 ’ 7 although the decrease is 
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much greater with the incorporation of radiation and two- 
phase flow losses 

Conclusions 

Based on the above analysis of an Al/hydrocaibon/0 2 
fueled rocket engine, the following results were obtained: 

1 . Exercise of the one-dimensional combustor model 
predicts that only moderate secondary atomization 
((3=5) is required to reduce overall propellant bum- 
out distance by 62% and final A1 2 0 3 residual di- 
ameter by 4 1%. 

2. Radiation losses for a 60 wt% gel, assuming |3=5, 
are estimated to be approximately 5 % of the energy 
released during combustion, resulting in a 1% de- 
crease in engine 1^. Secondary atomization may 
reduce radiation heat transfer losses, primarily due 
to decreases in the propellant burnout distance. For 
example, a fragmentation ratio of five results in a 
61 % decrease in radiation losses. Radiation losses 
were also found to be a function of gel composition 
and engine operating conditions, as a result of 
changes in chamber temperature. 

3. Two-dimensional, two-phase nozzle code results in- 
dicate that secondary atomization may have little ef- 
fect on nozzle two-phase flow losses. Specifically, a 
fragmentation ratio of 5 decreases two-phase flow 
losses by only 4 % compared to the above 62% re- 
duction in propellant burnout distance. Further- 
more, secondary atomization may have no effect on 
two-phase flow losses if particle coagulation, sur- 
face growth, and shear induced breakup are the 
dominant mechanisms controlling oxide particle 
size. 
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